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It is well known that the supersolid form of matter can arise in a single species of cold bosonic 
atoms in an optical lattice due to long range interactions. We present a scenario where a supersolid 
is induced in one of the components of a mixture of two species bosonic atoms where there is no 
long range interactions. We study a system of normal and hardcore boson mixture with only the 
former possessing long range interactions. We consider three cases: the first where the total density 
is commensurate and the other two where they are incommensurate . By suitable choices of the 
densities and the interaction strengths of the atoms, we predict that the charge density wave and 
the supersolid orders can be induced in the hardcore species as a result of the competing interatomic 
interactions. 
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The hallmark of the supersolid (SS) phase is the co- 
existence of the superfluid and the charge density wave 
(CDW); i.e. solid orders [l|,[2. This phase has not been 
observed unambiguously in experiments. However, in re- 
cent years several theoretical predictions of this phase 
have been made in different lattice systems 0, 0] ■ 
After the claim of the observation of this phase by Kim 
et al [1] which later could not be verified [13, [lH the 
search for the supersolid phase has become an active area 
of research. 

The pioneering observation of the superfluid (SF) to 
Mott insulator (MI) transition in an optical lattice using 
cold bosonic atoms which had been predicted by 

Jaksch et. al. [3] based on an earlier work by Fisher 
et. al. [3], has opened up new directions in the field 
of ultracold atoms. The possibility of tuning the inter- 
atomic interactions in optical lattice makes this system 
ideal for obtaining exoticphases of ultra cold atoms with 
long range interactions [3, [l^ . The observation of BEC 
in ^^Cr atoms which have large magnetic dipole moments 
in a trap [l^ in combination with the advancing research 
in optical lattice systems raises the hope of the observa- 
tion of the supersolid phase in the future. 

Theoretical studies of the possible existence of super- 
solids in mixtures of bosonic atoms and Bose- Fermi mix- 
tures have attracted much attention 17, mili. Mix- 
tures of ultracold atoms are very interesting because of 
the various competing interactions between the atoms 
leading to many exotic phases. It has been shown that 
for a single species ultracold atoms in a one dimensional 
lattice, the supersolid exists for large onsite as well as 
nearest neighbour interactions when the density of the 
system is not commensurate to the lattice and also larger 
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than half filling [3, In this Letter, we have con- 

sidered a mixture of two species bosonic atoms with 
one species consisting of normal and the other hardcore 
bosonic atoms. For the latter species, a single lattice site 
can be occupied by no more than one atom. This mix- 
ture can therefore be considered equivalent to a system of 
Bose-Fermi(spinless) atoms. We assume that the normal 
bosonic species exhibits long range interactions, but the 
inter atomic interactions in the hardcore species arc lim- 
ited only to onsite interactions. The model Hamiltonian 
for such a system can be written as: 
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Here and 6^, respectively, are the bosonic annihilation 
operator for atoms of a (normal) and h (hardcore) bosons 
localized on site i, nf = alai and = blbi represent its 
number operators and and are the hopping ampli- 
tudes between the nearest neighbours (ij). U"" {V"") are 
the on-site (nearest neighbour) intra-species repulsive in- 
teractions for the normal atoms. The inter-species ( be- 
tween normal and hardcore bosons) interaction is given 
by U""^. The hopping amplitudes (i", t^) and interac- 
tion parameters (C/°, V"", U"''') are related to depth of 
the optical potential, recoil energy and the scattering 
lengths [il,l2l|. The ratio U^'^/U'' as weU as U^'/V' 
can be varied over a wide range of values experimentally 
2^ . 2^ . In this work we consider t°- =t^ =t and set our 
energy scale by taking t = \. 

We identify various ground states phases of the model 
(H]), by calculating (i) the single particle excitation gap 
G1 for species a = a,b defined as the difference be- 
tween the energies needed to add and remove one atom of 
species a; i.e., Gf, = ELiNa + l,Nb)+ELiNa ~ l,Nb) - 
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2EL{Na, Nb) and (ii) the on-site number density defined 

by 



«) = ii'LN^N.lnfltlJLN^m)- 
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Here a, as mentioned before, is an index representing 
normal (a) or hardcore (b) bosons, with Na (Nb) corre- 
sponds to total number of a (b) bosons in the ground state 
li^LNaNb) of ^ system of length L with the ground state 
energy EL{Na, Nt,). The former is used to distinguish 
the gapless superfluid phase from the Mott insulator or 
the charge density wave phase, both having finite gap in 
their energy spectrum. In one dimension the appearance 
of the SF phase is indicated by for i oo. 

However, for a finite system is finite, and we must 
extrapolate to the L — > oo limit, which is best done by 
the finite size scaling of the gap 3, 111 ■ In the critical 
region 



Gt = L-'f{L/C) 
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where ^" is the correlation length for species a which di- 
verges in the SF phase. Thus plots of LGf^ versus the 
nearest neighbour interaction for different values of L co- 
alesce in the SF phase. On the other hand, when this 
trend does not follow, then the system is considered to 
be in the gapped, MI or CDW phase which is further dis- 
tinguished from each other via the CDW order parameter 
defined as 



O'^CDwiL) = \Y.^'^^^^nMK - p'^mLN.N,). (4) 



The existence of the solid order in the thermodynamic 
limit is verified from the finite value of Og^^y(L — > oo). 

We have employed the finite size density matrix 
renormalization group (FS-DMRG) method with open- 
boundary conditions to determine the ground state. This 
method has proved to be one of the most powerful tech- 
niques for studying ID systems 24, 25, 2^, 27 1. For the 
normal species, we have taken the maximum occupation 
per site as four (4). We allow up to 128 states in the 
density matrix of the left and right blocks in each itera- 
tion of the FS-DMRG calculations. The weights of the 
states neglected in the density matrix of the left and right 
blocks are less than 10~^. 

The charge density wave phase in bose systems is pos- 
sible when the density of bosons are commensurate with 
the underlying lattice. For example, the earlier stud- 
ies of the one-dimensional single species extended Bose- 
Hubbard model have shown the existence of the CDW 
phase for p = 1/2, 1 d H H llll. Later this study 
was extended in the case of two species extended Bose- 
Hubbard model, where the solid order is achieved for 
= = 1/2 by suitably varying the strengths of the 



nearest neighbour interactions [29[. Supersolid phase is 
then possible only moving away from these commensu- 
rate densities. In a recent study on a Bose-Fermi mixture 
with different hopping amplitudes, the supersolid phase 
has been predicted without the nearest neighbour inter- 
actions by doping the bosonic species [l3|. In order to 
achieve the supersolid phase in these systems, the inter 
atomic interactions have been carefully controlled such 
that the added bosons do not destroy the CDW phase by 
occupying sites that are already occupied. 

The recent study of a two species Bose mixture in a one 
dimensional lattice shows that phase separation occurs if 
the ratio U""^ /U"" is larger than unity [sO]- In order to 
avoid this condition, we consider U°' = U""^ = U and 
study the effect of on the ground state of model ([T|) 
for three possible combinations of densities; (i) p"' — p'^ — 
1/2, (u) p" = 1/2 and = 1/4 and (iu)p° = 3/4 and 

~ 1/2. In the first case, the total density of the system 
is commensurate, but in other two cases it is not. In all 
the above three cases we have taken J7 = 6, which is 
very large compared to the nearest neighbour tunneling 
amplitude t = 1. 

p° = p'' — 1/2 :- In this case the total density of bosons 
p = p°- + p^ = 1. For y = 0, the system is in the MI 
phase because the onsite intra species interactions, U"" = 
6, J7'' = oo and the inter species interaction, If^^ = 6 
are all greater than Uc ~ 3.4, the critical strength of 
the on-site interaction for the SF-MI transition in the 
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one-dimensional Bose-Hubbard model 
system continues to remain gapped as V increases. The 
gap corresponding to lattice size L for species a, G^ is 
plotted for different values of in Fig[TJ This figure 
clearly shows that G1^^ is finite for all the values of V"' 
that we have considered. However, the gapped phase at 
higher V is not a MI but a CDW, since 0^jjy^{L oo) 
is finite for V°- > ^ 1.0. In Figl2] we have plotted 
^CDwi^) versus 1/L for different values of V"". The or- 
der parameter Og.^^y(L) goes to zero for small values of 
and branches out for higher values indicating the on- 
set of the CDW phase. Op^^y(L) also exhibits a similar 
behaviour indicating that both the normal and hard-core 
bosons undergo a MI to CDW transition. The depen- 
dence of this transition on for the normal bosons is 
expected on the basis of an earlier work [i^. However, 
it was not obvious that the hard-core bosons would also 
undergo a similar MI to CDW transition as they lack 
long range interactions to exhibit density oscillations. 

The physical scenario when = p** = 1/2 is the fol- 
lowing: when y = 0, due to the strong repulsion be- 
tween the bosons, the system is in the MI phase with 
both normal and hard-core bosons uniformly distributed 
through out the lattice giving the average density of the 
total number of bosons at every site equal to one. As 
y increases, there is competition between the interac- 
tions, ?7°^, and and hence an atom of species a 
(normal bosons) cannot occupy the sites next to another 
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FIG. 1: (Color Online)Single particle energy gap G1 for nor- 
mal bosons as a function of 1/L for different showing finite 
gap in the thermodynamic limit. 
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FIG. 2: (Color online) Ocdw{L) versus 1/L for different V 
showing the onset of CDW order for V > VS ~ 1.0. 

atom of the same species, thereby forming a CDW phase. 
In addition, the atoms of the hardcore species b cannot 
occupy a site where there is either a hardcore boson or 
a normal boson because of the strong repulsive onsite 
interaction, t/"''. These physical conditions give rise to 
the intermingled CDW phase where the nearest neigh- 
bour sites are occupied by atoms of different species as 
shown in (Figl3]). It is interesting to note that the pres- 
ence of is sufficient to induce the solid order in the 
hardcore species in spite of the absence of any long range 
interaction between them. This type of induction of the 
solid order makes the other combinations of densities pre- 
sented below very interesting. 

p°- = 1/2 and p'' = 1/4 :- In this case the total bo- 
son density p — 3/4, and it is not commensurate with 
the lattice. In the present problem, we have not consid- 
ered long range interactions beyond the nearest neigh- 
bour, the commensurate densities are therefore intergers 
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FIG. 3: (Color online) The onsite number density for a and b 
bosons are plotted against the site index i taking V" = 5 for 
= 1/2, = 1/2. The black curve corresponds to < > 
and the red corresponds to < >. The node at the center 
is due to open boundary condition. 

or half integers. In a normal two species bosonic mix- 
ture with incommensurate density (e.g.p = 3/4), there is 
no transition from a SF to a gapped phase [30|. How- 
ever, such a transition does occur in a normal - hardcore 
boson mixture described by the model ([1]). The finite 
size scaling of the gap LG^ that is obtained from our 
FS-DMRG calculation shows, a transition from the gap- 
less SF phase to a gapped phase for the normal bosonic 
species in Fig|4l The critical value of is Vq ~ 3.0. 
However, the hardcore species remains in the SF phase 
showing no gap in the excitation spectrum. The calcu- 
lation of the CDW order parameters given in Fig. O for 
both the normal and the hardcore bosons show a finite 
OcDW in the hmit L ^ oo for > - 3.0. Thus 
the gapped phase of the normal bosonic species is iden- 
tified to be the CDW phase while the gapless phase of 
hardcore bosons is a supersolid since both sperfluid and 
CDW coexist. The CDW oscillations are similar to the 
case of Pa = Pfc = 1/2 as given in FigO 

= 3/4 and p'' = 1/2 :- The phase transitions that 
we have obtained for this case for the two species are 
in the reverse order as that of case (ii). The normal 
species shows a transition from the SF to the supersolid 
phase not showing any gap in the excitation spectrum 
(see FiglS]) and a finite CDW order parameter at ~ 
1.2 (see Figl?]). However, the hardcore species makes a 
transition from the SF to the CDW phase at the same 
critical point ^ 1.2. 

In recent experiments on dipolar atoms, it has been 
shown that the the ratio between the on-site interaction 
and the nearest neighbour interaction can be controlled 
by Feshbach resonance [1^. In the present work, we take 
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FIG. 4; (Color online)Finite size scaling of gap LGf^ is plotted 
as a function of V^" for different system sizes for = 1/2 
and p*" — 1/4. (a)The coalescence of the different curves for 
V° ^5 3.0 for the normal species shows the transition from the 
gapless SF phase to a gapped phase, (b) The coalescence of 
the different curves for all the values of V" for the hardcore 
species indicates the existence of the SF phase. 
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FIG. 5: OcDw{L) versus l/L for (a) normal and (b) hardcore 
species for V" ranging from 0.2 to 5.0 in steps of 0.4 is shown 
for = 1/2 and p'' = 1/4 . It is clear from the scaling that 
the O^DW becomes finite for ~ 3.0 for both the species 
indicating the transition to the CDW phase. 



FIG. 6; (Color online) Scaling of gap LGl is plotted as a func- 
tion of V" for different system sizes for = 3/4 and p'' = 1/2. 
(a)The coalescence of the different curves for all the values 
of V for the normal species shows that the system remains 
in the gapless SF phase. (Inset) Shows the G1 plotted against 
l/L for different V" ranging from 1 to 5 in steps of 1 which 
shows that the gap vanishes in the thermodynamic limit as 
the curves extrapolate to zero as L — > oo(b) The coalescence 
of the different curves for ^ 1.2 for the hardcore species 
indicates the transition from the SF to the CDW phase. 



Conclusion:- We have considered a system of normal 
and hardcore bosonic mixture with the normal species 
possessing long range interactions. By taking three dif- 
ferent sets of densities of both the species, we have in- 
vestigated the conditions that give rise to the supersolid 
phase in either or both the species. The main findings of 
this work is that by suitably tuning the nearest neighbour 
interaction strength V^, the solid order can be stabilized 
in the normal bosonic species and it can also be induced 
in the hardcore species as a result of the competition be- 
tween the J7", 11°''' and V°. This induction of the sohd 
order can lead to both the species being in the CDW and 
more interestingly one species in the CDW phase and the 
other in the supersolid phase depending on the choice of 
densities. By keeping the onsite repulsion for the normal 
species = 6 and varying the nearest neighbour in- 
teraction strength of the hardcore species V"", we obtain 
different interesting quantum phases which are listed in 
the table shown below. 



a fixed on-site interaction and vary the nearest neighbour 
interaction, which is equivalent to the experimental sit- 
uation where the ratio is varied. Our results suggest a 
way to obtain the supersolid phase, where we have re- 
stricted the control parameter only to the s-wave scat- 
tering length of one species which induces supersolidity 
in the other species. 



TABLE I: Different phases 
Pa pb softcore hardocre 

1/2 1/2 MI-CDW MI-CDW 
1/2 1/4 SF-CDW SF-SS 
3/4 1/2 SF-SS SF-CDW 
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FIG. 7: Finite size scaling of Ocdw for (a) softcore and (b) 
hardcore species for V"' ranging from 0.8 to 2.4 in steps of 0.2 
is shown for p" = 3/4 and = 1/2 . It is clear from the 
scaling that the Ocdw becomes finite for V ~ 1.2 for both 
the species indicating the transition to CDW phase. 



If the individual density of a species is commensu- 
rate then there exists a SF-CDW or MI-CDW transition 
where as if it is incommensurate then the transition is 
from SF to SS. It is interesting to note that for incom- 
mensurate fillings the transition from the SF phase to the 
CDW or SS phases occur at the same critical value Vq 
for both the species. 
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